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We study the Wheeler-DeWitt quantum cosmology of a spatially flat Friedmann cosmological 
model with a massless free scalar field. We compare the consistent histories approach with the 
de Broglie-Bohm theory when applied to this simple model under two different quantization schemes: 
the Schrodinger-like quantization, which essentially takes the square-root of the resulting Klein- 
Gordon equation through the restriction to positive frequencies and their associated Newton- Wigner 
states, or the induced Klein-Gordon quantization, that allows both positive and negative frequencies 
together. We show that the consistent histories approach can give a precise answer to the question 
concerning the existence of a quantum bounce if and only if one takes the single frequency approach 
and within a single family of histories, namely, a family containing histories concerning properties of 
the quantum system at only two specific moments of time: the infinity past and the infinity future. 
In that case, as shown by Craig and Singh [TJ, there is no quantum bounce. In any other situation, 
the question concerning the existence of a quantum bounce has no meaning in the consistent histories 
approach. On the contrary, we show that if one considers the de Broglie-Bohm theory, there are 
always states where quantum bounces occur in both quantization schemes. Hence the assertion that 
the Wheeler-DeWitt quantization does not solve the singularity problem in cosmology is not precise. 
To address this question, one must specify not only the quantum interpretation adopted but also 
the quantization scheme chosen. 



I. INTRODUCTION 

It has been claimed in many papers (some few exam- 
ples are Refs. pU!]) that the Wheeler-DeWitt approach 
to quantum cosmology [6, 7 does not solve the singular- 
ity problem of classical cosmology. This quite strong and 
general assertion can be criticized in many ways. First 
of all, this claim is usually based on calculations on a 
very simple model, namely, a free massless scalar field 
in Friedmann models, which is of course a very narrow 
subset of cosmological models. Secondly, the quantiza- 
tion program which is carried out on those papers is very 
particular: the Wheeler-DeWitt equation of these mod- 
els are Klein-Gordon like equation, and the procedure 
is to extract a square root of it and work in a single fre- 
quency sector. Note that this is not mandatory, and there 
are other ways to deal with the Klein-Gordon equation 
working in the two frequency sectors with a well defined 
inner product, as it can be seen in Refs. |9]. Fur- 
thermore, inside this single frequency sectors, the analy- 
sis is a priori restricted to quantum states which are ei- 
ther right-moving (expanding classical solutions) or left- 
moving (contracting classical solutions) [2H5] . Finally, in 
some cases, the interpretation of the quantum states is 
not described clearly [3H5] . Expectation values of the vol- 
ume operator are calculated, but what do they mean? Of 
course they are not averages of external observers mea- 
surements of the volume operator as long as we are deal- 
ing with cosmology. Hence, are these expectation values 
interpreted along the lines of the many worlds interpreta- 
tion [10], the consistent histories approach [TT], or some- 
thing else? Among the above mentioned papers, only 



Refs. [TJ [2] identify which interpretation is being used. 

In the framework of a single frequency quantization of 
the simple model presented above, reference [JJ presents 
the most rigorous and precise approach to the question 
of the existence of a quantum bounce in the Wheeler- 
DeWitt quantization. The interpretation they adopt is 
precisely defined (the consistent histories approach), and 
the conclusion is that, if one takes family of histories with 
properties defined in just two moments of time, the infin- 
ity past and the infinity future, then the probability of a 
quantum bounce is null for any quantum state, including 
superpositions of right and left-moving states. This is a 
remarkable result. 

The aim of this paper is to discuss the results of 
Ref. pQ with care, and contextualize it in the framework 
of other interpretations of quantum mechanics, namely, 
the de Broglie-Bohm theory, and other quantization tech- 
niques, as the two frequencies (Klein-Gordon) approach 
of Ref. 0. 

We will first show that, in the single frequency ap- 
proach using the consistent histories interpretation, fam- 
ilies of histories containing properties defined in one or 
more moments of time, besides properties defined in the 
infinity past and in the infinity future, are no longer con- 
sistent, unless one takes semi-classical states, which of 
course corresponds to histories without a bounce. This 
means that in the framework of these families of histories 
one cannot answer whether quantum bounces take place 
because histories involving any genuine quantum states 
are inconsistent. Hence the consistent histories approach 
is silent about quantum bounces happening in family of 
histories with more than two moments of time. Further- 
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more, we will show that in the induced Klein-Gordon ap- 
proach there are no consistent family of histories involv- 
ing genuine quantum states. Again, the question about 
the existence of quantum bounces has no meaning in the 
induced Klein-Gordon approach. 

On the contrary, if one considers the de Broglie-Bohm 
theory, where trajectories in configuration space are con- 
sidered to be objectively real (the so called Bohmian tra- 
jectories), there is a vast number of examples of non- 
singular models. In this case, one can show that in the 
two quantization procedures mentioned above, there ex- 
ist plenty of bouncing trajectories which tend to the clas- 
sical cosmological trajectories when the volume of the 
universe is big. Similar results have already been ob- 
tained in earlier works (see Refs. [HJ-JTB]). Hence, the 
existence of quantum bounces in the Wheeler-DeWitt ap- 
proach depends strongly on the quantum interpretation 
one is adopting, and on the quantization procedure one 
is taking. 

The paper will be divided as follows: in section II we 
will present the minisuperspace model we will work on, 
and summarize the results of Ref. [T]. In section III, 
still in the framework of Ref. |T| , we will consider fami- 
lies of histories with properties defined in three or more 
moments of time, and we will show that they are not 
consistent for general quantum states. In section IV we 
will summarize the de Broglie-Bohm theory applied to 
quantum cosmology and we will present an infinitude 
of Bohmian trajectories which are non-singular and ap- 
proach the classical limit for large volumes of space. In 
section V we show that there are no families of consis- 
tent histories in the Klein-Gordon quantization for gen- 
eral quantum states. However, also in this quantization 
scheme, there are many Bohmian trajectories which are 
non-singular. We end up with the conclusions in section 
VI. 



to have a unit volume. Denoting a :— log a, we rewrite 
the Hamiltonian in a more convenient form 



,-3a 



H = 



AttG 



-Pa 



0. 



(3) 



Note that our choice of variables is different of the one 
used in pQ. The conclusions are nevertheless unchanged. 
The pair (a.p a ) is still canonical. After quantization 
the momenta promoted to derivative operators and the 
classical constraint H sa becomes the Wheeler-deWitt 
equation 



(4) 



defined on the kinematical Hilbert space L 2 (M 2 , dadcf>) . 
We see that the Wheeler-deWitt quantization of this sim- 
ple model is equivalent to the quantization of the Klein- 
Gordon equation in Minkowski spacetime. The standard 
procedure is to separate the positive and negative fre- 
quency modes and quantize them independently. Taking 
the square-root of the constraint, we get 



±id <t> ^{a.4>) = V0 *(a,< 



with 







4vrG 



(5) 



(6) 



The action of v© is best seen on Fourier space. Con- 
sider the set of eigenf unctions 



e k (a) = (a\k) - 
such that e,t = to 2 ek, with 



>2it 



(7) 



II. THE CRAIG AND SINGH RESULT 

In this section we shall briefly develop the Whceler- 
DeWitt quantization of a FLRW model with flat spa- 
tial sections following closely reference [1]. We assume 
that the matter content of the universe can be described 
by a massless scalar field and we consider the consistent 
histories approach of [IT] to explore the existence of a 
quantum bounce. The Hamiltonian constraint of this 
minisuperspace system reads 



H 



3 a 2a 3 



0. 



(1) 



where a is the scale factor, (f> the scalar field and p a ,'j 
their conjugate momenta, given by 



Pa = 



4ttG 



aa , Pcf, = a 



(2) 



On the derivation of the Hamiltonian above the lapse 
function was set to one and the fiducial cell considered 



k\. 



The physical scalar product is given by 



($|*) := I da# 

J 4>=<f> 



(8) 



(9) 



and it is independent of the time 4>q on which it is de- 
fined. Positive and negative frequency sectors are orthog- 
onal with respect to this scalar product. Restricting to 
the positive frequency sector, evolution is given by the 
propagator 



C/(0-0 o ) = e^ 



(10) 



There are two Dirac observables. The first one is p^, 
which is an invariant of the model. The second one is 
a relational observable [IT]- For any operator A, which 
does not need to commute with the Hamiltonian, one 
constructs the corresponding relational observable 



lU*(<£)) = [^-0o)i|*(0o)>, 



(11) 
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which implies 



M<h =u(c/> -cf>yAu(<t>o-4>)- 



(12) 



Applying this to the operator a, the action of the Dirac 
observables on the positive frequency sector is given by 



a| 0o *(a, <f>) = U{(j) - o )a*(a, fa) 
p^(a,(j)) = h\f@y{a,4>). 



(13) 



Note that the action of the Dirac observables preserves 
the positive and negative subspaces, which is consistent 
with the approach taken. 

With the physical Hilbert space specified, we need to 
construct now the set of histories and the corresponding 
decoherence functional. 

Following Hartle's approach [11] we are interested in 
defining a decoherence functional for a set of histories. 
The decoherence functional is defined as 



(14) 



where the branch wave function is given by 

**:=C£|*). (15) 

In the above formula, \l/ is a given initial state and Ch 
is the class operator defining the history h, given by a 
product of projectors 



C h := F 



(*1 



(in). 



(16) 



P^x {t) projects onto the subspace for which the kth 
eingenvalue of the observable O at time t takes values in 
the interval AA&. Here we are using Hcisenberg operators 
for the projectors 



(17) 



For the case at hand we consider the observable given 
by the scale factor a, or a, with relational time <j), and 
we will denote projectors simply by P&a i {4'i)- The time 
independent projector is given explicitly by 



PAa = I da\a)(a\ 

I Act 



(18) 



where the ket \a) is defined in ^ and the normalization 
is such as to make this basis orthonormal. Let us further 
compute the expectation value of the evolution operator 
in this basis, as it will be useful in the next section 



(at\U(4>)\a) 



dk 
27 



where we used the resolution of the identity on Fourier 
space 



dk\k)(k\ = 1. 



(20) 



The set of histories considered in [Tj are composed by 
two times, corresponding the past infinity (0 — > — oo) and 
future infinity (4> — > +oo). The histories are separated 
in those where a is bigger or smaller than a given fixed 
fiducial value a*. There are four possible histories ac- 
cording to these possibilities, described by the following 
class operators 

Cs-s(-°0,°o) = PAa 1 (-Oo)P A Q 2 (oo) 

C S -b(-oo,oo) = P Aqi (-oo)P Aq2 (oo) 
C B -s(-oo,oo) = Pa Qi (- 00 )- p Aq 2 (oo) 

Cb-b{~00,00) = PA ai (~°°) P Aa 2 (oo) 

where S and B subscripts denote, respectively, the do- 
mains of the scale factor arbitrarily close to the singular- 
ity or arbitrarily big. 

Let us now check that this set of histories is consistent. 
Noting that PaoPaq = 0, the only non-trivial terms are 
d{hs-B,hs-B) and d(hs-s,hB-s)- Consider for exam- 
ple the first of these terms 

d(h S -B,h B -B) = (*I^A Ql (0l)^Aa,(02)PAa 1 (^l)|*). 

(21) 

Let us as a first step study the behavior of Paq(0)| x I') 
and P/± a (4>) for tj> ~ > ±oo. We borrow the results without 
proof from pQ. We have that 

lim P A a(4>)\y) = 

q>— >-+oo 



lim P Aq (0)|*) 
lim P Aa (<t>)\V) 

>— y-\-oo 



\*r) 



lim P Aq (0)W = |* l) 

q)—^ — oo 



(22) 



In the equation above we have used the left/right- 
moving decomposition of the wave function 



1 



'2tt 



dk^{k)e lka e lUHt} 



&*Ofc)e lfe(Q - 0) 



dk^(k)e 



(23) 



where v r := a — (ft, vi := a + <p and we dropped the factor 
of y/AnG/3 in u. 

Since right and left moving sectors are orthogonal, the 
term in (21 ) is zero, as is the other term, and the deco- 



herence functional is diagonal for this set of histories. 

Craig and Singh go on and restrict the histories to 
those giving rise to a bounce and those that have a sin- 
gularity and prove that the probability for a bounce is 
zero. 

As noted by the authors in pQ, the decoherence func- 
tional in the case with more than two times is not neces- 
sarily diagonal, unless the wave functions are taken to be 
semiclassical. In the next section we will show explicitly 
that already in the case of three times the decoherence 
functional is not diagonal. 
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III. HISTORIES WITH THREE TIMES 

Let us now see whether we can obtain a family of con- 
sistent histories when we ask about properties concerning 
the size of the universe in a third moment of time between 
<t>\ — > — oo and 2 —> +00. Thus we want to address the 
question whether in an arbritary intermediary cj> time the 
scale factor of the universe is in the interval (—00, a*), 
or on its complement (a*, 00). 

The new family has now eight histories associated with 
the following class operators 



C S - 


-Aa- 


-s( 


<t>\A 


,<h) 


= PAocA 


4>l)PAa((f>)PAa 2 { 


02), 


C s - 


-Aa- 


-s( 


faA 


,02) 


= PAccA 


0l)^Aa(0)^W 


02), 


C s - 


-Aa- 


-b{ 


fa,<t> 


,02) 


= PAccA 


0l)PA a (0)PA Q2 ( 


02), 


C s - 


-Aa- 


-b( 


i>\A 


,02) 


= PAccA 


0l)^A a (0)^A Q2 ( 


02), 


C B - 


-Aa- 


-s( 


i>\A 


,02) 


= Paoh ( 


0l)^Aa(0)iW 


02), 


C B - 


-Aa- 


-s( 


i>i,<t> 


,02) 


= PAaA 


0i)i 5 A a (0)^W 


02), 


C B - 


-Aa- 


-b{ 


i>\A 


,02) 


= PaoA 


0i)Pa q (0)Pa Q2 ( 


02), 


C B - 


-Aa- 


-b( 


i>i,<t> 


,02) 


= PAaA 


^PAai^P-AaA 


02124) 



with S and B having the same meaning as before being 
close to the singularity or arbitrarily big. 



Each of these class operators is associated with a 
particular history. For instance, the class operator 
Cs-Aa-s(0i, 0, 02) is associated with the history where 
the universe was singular at (f>i — » —00, has a size in a 
domain Aa at the finite time 0, and it will be infinitely 
large at 02 — » 00. 



One must now see whether this new family with eight 
histories is consistent or not. As we have seen above, 
one must calculate the decoherence functional d(h, h') for 
the histories associated with the class operators shown in 



Eq. (24) 



It is easy to show that, in general, d(h, h!) is not ap- 
proximately zero. For that, let us calculate the decoher- 
ence functional for the histories associated with the class 
operators Cs_ac-b(0i, 0, 02) and C s _ Aa _ B (<j)i, 0, 2 ). 
In this case the functional reads 



d(h s _ Aa _ B ,h S -Aa-B) = Tr(P AQ2 (^ 2 )P AQ (0)P AQ1 (^ 1 )|*)(*|P Aai (0 1 )P Aa (0)P AQ2 (0 2 )), 

= (^|PAa(0)PA Q2 ^A Q (0)l*fl), 



(25) 



where we have used that P Aai (0i — > — oo)|\I r ) = \^n)- Then, dropping constant factors in front of the integrals and 
using eq. (19), we have 



/a* />OG />OG 

da"*>;') / da*(v r ) da' < a"|[/(0 - 2 )|a' >< a'|f/(0 2 - 0)|a > 
-00 Ja* Ja 2 * 



lim / da"9*(v'^) da^{v r ) da' ( dk' e -W + / dk' e -**V-""+*-*») 



dk e -iH<*-<*'+<i>-<h) + dk e -*H<*-<*'+4>2 



(26) 



We shall analyse each of its four terms separately. The term formed by the product of the first with the third can 
be written as 

lim / e^-^dvU *(v r )da da"**«) / dk' e tk ' v ' dk e~ lkvi = , (27) 

while the product of the second with the fourth gives 



lim 



a 2 * — (p2 
oo pa* 



*<(*-*') / V( Vr )da / da"y*{v") / dk' e lk '< / dk 



-ikv r 



ip.v. 



V(v r )da I da"^*(K) I dk e- ik{ <- v ^ 
„/¥(v r )**«) 



$t(v r )dv r I d<**«) 

—00 

a* — <f> 



dv r 



dv'l 



(28) 
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where in the last equality we have used the fact that the integral in v r and v" are in disjoint domains, and hence the 
part involving the delta function is null. For the definition of the principal value, noted above as p. v., see appendix 
[A") The sum of the cross terms involving the products of the first with fourth terms, and the second and third reads, 
respectively, 



lim 



*(v r )da I da"¥*«) / dv'J dk' e - ife 'K-<+ 2 (^-«] / dke~ ik{ <- v ^ 

> Ja 2 *-<f>2 \Jo JO 



/>oo poo \ 

+ j dk e -**[«r-»;-a(* a -«] j dk > e -ifc'«'-oj _ 



(29) 



Let us concentrate on the first term. The other term follows the same reasoning. Performing the integrals on k and 
k' , one gets 



lim / ty(v r )dv 



d<**«) / dv' r 

1 Jct2* — (f>2 

7rtf«-< + 2(^,-0))-ip.v/— ' 



vL - v" + 2 



The typical terms are proportional to 



lim 

j>2-HX> 



lim p.v. 



2(02-0)) 



a* — <f> poo 



x 
a* — 



a* — <f 
oo 



*(fr) 



7T(5(u' — u r ) — ip.v. — 

\ nv 



V" — V r 



dv"V*(v") / dv r V(v r ) / da' 

-oo Ja* — (b J a.2* 



(a' - «» + 2 - 20)(a' - w r - 2 )y ' 



(30) 



(31) 



The terms with &*(v" — 2(02 — 0)) are zero in the limit 02 — > oo because ^ is assumed to be square integrable, while 
the terms with a 2 in denominators are obviously null in this limit. 



The only non-null term of Eq. ( |26| ) is Eq. (28). Hence, the final result for this off-diagonal term of the decoherence 
functional is 



d(h 



S-Aa-Bi n S-Aa-B 



) OC — ip.v. 



dv r 



a* — (fi 



dv'J 



*(tv)**K.') 



(32) 



which is not null in general. Due to the disjoint domains of integration, this result can be approximately zero if 
and only if ^(v r ) is concentrated around some fixed value of v r . The classical trajectories are given by v r — const 
or vi — const. Therefore, a wave function sharply concentrated around some fixed value of v r must describe a 
semiclassical state. It is straightforward to show that other off-diagonal terms of the decoherence functional, e.g., 
dQiB-Aa-Si ^B-Aa-sli are approximately zero only if the wave function ^(v;) is concentrated around the other class 
of classical trajectories vi = const. 



Concluding, the family of histories described by the 
class operators ( 24 1 can be made consistent only for semi- 
classical states. In that case, of course, the probabil- 
ity of occurrence of a quantum bounce is null, as be- 
fore, but the reason for that comes from the fact that 
we are not allowed to calculate probabilities in a fam- 
ily of cosmological histories where quantum effects arc 
relevant. Probabilities are calculable only for semiclas- 
sical histories, where bounces cannot occur. More gen- 
erally, if quantum effects are important in any family of 
cosmological histories, under the consistent histories ap- 
proach one cannot ask any questions about properties of 



the Universe at an arbitrary finite 0. This is of course a 
limitation on the applicability of the consistent histories 
approach to cosmology, at least for the present simple 
model. We are simply prohibited to study the quantum 
properties of a cosmological model, unless one considers 
just two moments of its history, at ± oo, and noth- 
ing more than that. Are there any other approaches to 
quantum cosmology where one can go further? 
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IV. THE DE BROGLIE-BOHM THEORY 
APPLIED TO QUANTUM COSMOLOGY 

A quantum theory that can be consistently im- 
plemented in the quantum cosmology scenario is the 
de Broglie-Bohm quantum theory (see Ref.'s [H]-[23] for 
details). Considering minisuperspace models, which have 
a finite number of degrees of freedom, the general form 
of the associated Wheeler-De Witt equation reads 



dq p dq a 



UfaMq) = 



(33) 



where f p(y {qp) is the minisuperspace DeWitt metric of the 
model, whose inverse is denoted by f p<7 (q P )- By writing 
the wave function in its polar form, \1/ = R e lS , the com- 
plex equation ( 33 ) decouples in two real equations 



1 



dS dS 



2 f P M- d(lp dqo 



+ U(q li )+Q(q lt ) = 



d 



R z 



dS 



dq p \ dq a 



where 



= 



d 2 R 



(34) 



(35) 



(36) 



is called the quantum potential. The de Broglie-Bohm in- 
terpretation applied to Quantum Cosmology states that 
the trajectories q u ( t) ar e real, independently of any obser- 
vations. Equation (f34J) represents their Hamilton- Jacobi 
equation, which is the classical one added with a quan- 
tum potential term Eq. ( 36 ) responsible for the quantum 
effects. This suggests to define 



V 



as 

dq P ' 



(37) 



where the momenta are related to the velocities in the 
usual way 



1 N dt 



(38) 



In order to obtain the quantum trajectories, we have 
to solve the following system of first order differential 
equations, called the guidance relations 



dS(q P ) 1 
dq p 1 N q °- 



(39) 



The above equations ( 39 ) are invariant under time 



reparametrization. Therefore, even at the quantum level, 
different time gauge choices of N(t) yield the same space- 
time geometry for a given non-classical solution q a (t) . In- 
deed, there is no problem of time in the de Broglie-Bohm 



interpretation for minisuperspace quantum cosmological 
models [25] . However, this is no longer true when one 
considers the full superspace (see [351 [25]). Notwith- 
standing, even with the problem of time in the super- 
space, the theory can be consistently formulated (see 

[SUES]). 

Let us then apply this interpretation to our minisuper- 
space model. The Wheeler-DeWitt equation then reads 



a 2 * a 2 * 

da 2 d(f> 2 







(40) 



Comparing Eq. ( 40 1 with Eq. ( 33 ) , we obtain, from 
Eqs. p4| and p5l, 



as 
da 



dS 

a/> 



+ Q{q„) = 



aa \ oa 



where the quantum potential reads 

* d 2 R _ cPR 

da 2 dd> 2 



Q{a,(j>) 



1 
R 



The guidance relations ( 39 1 are 



dS 
do 



e 3a a 
N 



as r 

~d4> 



3a, 



N 



(41) 



(42) 



(43) 



(44) 



(45) 



We can write equation Eqs. (41) in null coordinates, 
1 



vi := 



V2 



(a + 0) 



1 / 

a := -j= (vi+v r ) 



1 



V2 



(vi - v r ) (46) 



yielding, 



d 2 



*(uj,tv) = 



dvidv r , 
The general solution is 



(47) 



(48) 



where F and G are arbitrary functions. Using a separa- 
tion of variable method, one can write these solutions as 
Fourier transforms given by 



*(«J,«r) 



dkU(k) e lkv ' + I dkV(k) e***- 



(49) 
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with U and V also being two arbitrary functions. If one 
restricts the wave function Eq. ( 49 ) to left or right moving 



components only, the quantum potential will necessarily 
be a function of either vi or v r , and hence it will be null 
(see Eq. (43)). In this case, only classical trajectories, 
which arc of course singular, are allowed. Hence, avoid- 
ance of singularities is possible if and only if the wave 
function Eq. ( 49 1 depends on both left and right moving 
components. 



Under restriction to positive frequency solutions, one 



gets a subclass of the general solution Eq. ( 49 ) 



iS(vi,v r )= dk^(k) e ikv ' + dk^(k) e lkv 



(50) 

From the guidance equations (44) and (45), one obtains 
that 



da 



dS/da 
dS/d<t> 

dk \ dk! 
ii Jo 



*(fc)**(fc') e lVl{k - k ' ] - e ~*«r(k-k!) 



(k + k') 



iv r k —ivik' 



¥(fc)**(-fc') e m ' k ] e lVrk ' 



(k-k')\/ 



oo /*oo 

dk / dk' 
(i Jo 



+ *(-fe)**(fc') e- m " k e~ lVlk ' +^(k)^*{-k') e lVlk \ e™" fe ' 



(k + k') 



(51) 



Let us analize Eq. |5lJ) in the limits v r — > ±oo or vi — » ±oo. When i? r — > ±oo, the integrals involving 
J °° dkty(k) e mrk , J °° dkk^{k) e lVrk and similar correspond to a Fourier transform of square integrable functions 
which are null when evaluated at v r — > ±oo. Hence we obtain from Eq. (51 1, in this limit, that 



da 



= -1 



+ tfi = vi = const 



(52) 



For vi — > ±oo, an analogous reasoning yields 



da 



a — 4> 



const 



(53) 



Hence, in the regions v r — > ±oo and u; — > ±oo, the Bohmian trajectories emerging from Eq. (51) are the classical 
trajectories irrespectively of the wave function. 

We will now see, however, that there are a huge class of states where the Bohmian trajectories are not classical in 
other regions of the (a, cf>) plane. For instance, when ^(k) is even on k, Eq. (51 ) reads 



da _ J °° dkf™ dk'^{k)y*(k') e^( fc -' £ '){sin[a(fc-fc')](fc±fc , )+sin[a(fc + /fc')](fc-^)} 
# ~ ~ l J °° dk J °° dfc'*(fe)**(fc') e^( fe - fe '){ C os[a(£; - Wj] + cos[a(fc + k')}}(k + k') 



(54) 



Note that Eq. (54) is anti-symmetric under the change a — > —a and also da/d<j) — at a — 0. Consequently, 



the Bohmian trajectories that start at v r — > oo (infinitely big universe) cannot cross the line a = and goes to the 
singularity at v r — > —oo in the same way as the classical trajectories do. These Bohmian trajectories are non-singular. 
On the other hand, if they start at the singularity in vi —> — oo, they cannot become infinitely big at vi — > oo. 

Note that this result is in opposition to the consistent histories conclusion. These Bohmian trajectories describe 
exactly what the consistent histories approach has claimed to be impossible, namely, universe histories that start 
infinitely big in the far past and go infinitely big also in the far future. In fact, the even states shown above within the 
de Broglie-Bohm scenario violate the consistent histories description for all trajectories. As just argued, there is no 
single trajectory that can start infinitely big in the far past and goes to a singularity in the far future or the reverse. 
Hence, it is certain that there exist non-singular Bohmian trajectories. 



One can also obtain bounces in the situation where Vf^fc) is not only even on k but it is also real. In that case 

(55) 



Eq. ( 54 ) simplifies to 

da J °° dk J °° dk'^(k)^*(k') sin[0(fc - k')]{sm[a(k - k')}{k + k') + sin[a(fc + k')](k - k')} 



d4> f™dk J °° dk'V(k)$* (k 1 ) cos[0(fc - fc')]{cos[a(fc - k')] + cos[a(fc + k')]}(k + k') 

where we have used the fact that only even integrands can survive. This can be seen by performing a coordinate 
transformation in k space, 

u := — = (fc + k 1 ) k '.— —= (u + w) 

w:=-j=(k-V) k':=-j=(u-w), (56) 

changing the integral domains accordingly, J °° du J*" dw, and noting that ^/(u + w)^(u — w) is even under the change 



w — > —w. Note that now Eq. (55) is anti-symmetric under the change <p — > —(f) and again we have that da/d<j) = 
but at <j) = 0. Hence, the Bohmian trajectories must certainly present a bounce when it crosses the line (f> = 0, and 
if they start at v r — > oo in a classical contraction from infinity, they must necessarily end at vi — > oo in classical 
expansion to infinity, realizing a bounce at <f> = and never reaching the singularity, in a symmetric trajectory in 
On the other hand, if they start at the singularity in u; — > — oo, they come back to the singularity at v r —> — oo, with 
the turning point taking place at (f> = 0. Again, contrary to the consistent history conclusion, any universe history as 
described by these Bohmian trajectories coming from infinity must go back to infinity, and any Bohmian trajectory 
coming from the singularity must go back to the singularity. 

Note that the line a = 0, where these non-classical behaviors are strong, corresponds, in our units, to a p ij ys = l p \, 
where is the Planck length. All these features can be seen numerically with a particular example. Let us take, for 
instance, 

*(fc) = e-a' 1 !-'*) 2 /^ (57) 

with a << 1 and d > 1. This is a real and even ^(fc), consisting of two sharply peaked gaussians centered at ±d. 
The wave function reads 

/•oo />0 

®(vi,v r ) = / dk^(k) e ikv ' + / dkV(k) e ikVr 

JO J -oo 

/•OO 

dk e -( fe - d ) 2 /- 2 e lkv < + / dk e -( fc + d ) 2 /- 2 e lkv - . (58) 



The Bohmian trajectories associated with this wave 
functions can be seen from figure 1. We can distinguish 
two kind of trajectories. The upper half of the figure 
contains trajectories describing bouncing universes while 
the lower half corresponds to universes that begins and 
ends in singular states ( "big bang - big crunch" universe) . 

In general, there is also the possibility of trajectories 
describing cyclic universes. In Ref. [12], it was consid- 
ered Bohmian trajectories associated with wave functions 
similar to the above one, but without the restriction to 
positive frequencies only. Considering both positive and 
negative frequencies, there are oscillatory trajectories in 
cf>. In this case, if one wishes to interpret <j) as time, 
this corresponds to creation and annihilation of expand- 
ing and contracting universes that exist for a very short 
duration. This fact suggests that one cannot understand 
Eq. ( 42 1 as a continuity equation for an ensemble of tra- 



jectories with a distribution of initial conditions given 



by R 2 in this Bohmian approach with guidance relations 
defined as in Eqs. (44) and (45). In fact, this interpre- 



tation of a continuity equation would be possible only if 
Eq. ( 42 ) could be reduced to the form 



dR 2 







dcf> da 



R 



, da 
~d4> 



0, 



with da/d(j> given by Eq. (51 ) 



da dS/da 
d~4> = ~ dS/d(f> ■ 



(59) 



(60) 



It can be shown, using Eqs. ( 41|42|44|45 ), that this is 
possible if and only if 



as d 2 S _ dS d 2 S 
da dad<p d 



(61) 
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FIG. 1: The field plot shows the family of trajectories for the 
planar system given by ( 44 1 ( 45 1 for the wave functional ( 58 1. 



Two of them that describe their general behavior are depicted 
in solid line: the first representing a bouncing universe while 
the second one corresponds to a universe which begins and 
ends in singular states ("big bang - big crunch" universe). 



we loose the notion of probability. 

It is worth remarking that, in addition to the usual de- 
bate related to different approaches to describe a quan- 
tum system, every quantum model of the universe has 
to face a non-trivial problem. Quantum cosmology deals 
with a single system, which forbids us to repeat experi- 
ments, hence posing peculiar issues associated with the 
physical meaning of any kind of probability in this con- 
text. Thus, the lack of probabilistic predictions in quan- 
tum cosmology should not be taken a priori as a defi- 
ciency of the formalism. On the contrary, one should 
carefully analyze if one can consistently extract informa- 
tion and predictions from the model without a notion of 
probability. 

We should stress that one can recover the probabilistic 
predictions in quantum cosmology using the de Broglic- 
Bohm theory when one implements a more complex mod- 
eling of the universe, adding new degrees of freedom 
(here, we have only one degree of freedom). In that case, 
a probability measure naturally appears in the quan- 
tum description of the sub-systems of the universe (see 
Ref. [30] for details), and the usual Born rule can be re- 
covered. In that case, for the sub-systems, the consistent 
histories and the de Broglie-Bohm approaches should co- 
incide. 



which implies that <p = 0, stating that is a monotonic 
function of coordinate time. This is a strong restrictive 
condition, which cannot be satisfied by general quantum 
states as the ones discussed above. In general, 0^0. 
Hence, Eq. (42 1 cannot be interpreted as a continuity 
equation in (f> time for the ensemble of trajectories given 
by Eq. (60) with distribution R 2 , even in the single fre- 



quency approach where one has a Schrodinger-like equa- 
tion. 

If, however, one insists in interpreting </> as the time 
variable, then one would have to face the situation of 
creation and annihilation of universes which is a typical 
feature of relativistic quantum theory. Accordingly, the 
lost of a continuity equation for R 2 can be associated 
with the non-conservation of the number of trajectories 
of this ensemble in the <f> time. 

Normally the de Broglie-Bohm theory of a 
Schrodinger-like equation furnishes, besides the quantum 
trajectories, a probabilistic measure for these trajecto- 
ries. This is not the case here since the kinetic term in 
the present Schrodinger-like equation is not canonical, 
hence, it is not of the form gij{x)p l p> , where gij{x) has 
an euclidean signature. 

Concluding this section, in the consistent histories ap- 
proach we may have the notion of probabilities but we 
are not allowed to investigate non-classical properties of 
the universe in any finite </> time, or to have more than 
two snapshots of any non-classical universe. 

On the contrary, in the de Broglie-Bohm theory we 
can investigate the entire evolution of the universe, but 



V. THE KLEIN-GORDON APPROACH 

In the Wheeler-De Witt equation for a free massless 
scalar field, one can define its square-root and construct 
a Schrodinger-like equation as we have discussed in sec- 
tions [TT] and |III| However, there are other quantization 
schemes where the restriction to a single frequency sec- 
tor is not necessary. A promising alternative approach to 
quantum cosmology using the consistent histories quan- 
tization is to consider the full Klein-Gordon equation. In 
this approach both energy sectors, positive and negative, 
are simultaneously taken into account but the Hilbert 
space is defined with a different inner product (see [5] 
and references therein). 

Following closely Ref. [5] , one can define the eigenstates 
associated to the position operator as 



\x) 



1 f°° dk 
1 f°° dk 
= + \x-) , 



,i\k\4>-ika\fo\ 

e _ *'l*l^ - ** a |jfc_ 



(62) 



where \k±) are eigenstates of the k operator such that 
k\k±) = k\k±) and k Q \k±) = ±\k\ \k ± ). Note that the 
position eigenstates are not orthogonal, i.e. 

(x\x'} = G^+\x,x') + G^(x,x') 
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where 

G (±) (x,x') 



1 

27 



dk 
2\k\ 



^i\k\ (0— (f)')±ik(a — a') 



(63) 

are respectively the positive and negative Wightman 
functions. The positive and negative position eigenstatcs 
satisfy a completeness relation that reads 



1 



da 



(\x + ) d <f,(x + \ - \x_)dj,(x_\) 



Given these position eigenstates we can define the in- 
duced Klein-Gordon inner product as 



i l da ['$*d dl <£>+ - %*_dl$- 



(64) 



where \l/± (a, 0) denotes the positive (negative) frequency 
solutions of the Klein-Gordon equation which are given 
by the projection of the wave function in the position 
eigenstates. Recalling that vi = a + and v r = a — 0, 
we have 



<F+(0, a) = (x+\V) 
1 



/ dk e lkv ^+(k) + \ 

Jo J-i 



V2tt 

--*Uv r ) + ¥ + ( Vl ) 



dk 



ikvi 



(65) 



and 

^-(0, a) = 
1 



'2tt 



(i 



dk e lkv 'V-(k)+ / dk e lkVr V-(k) 



with 



*±(fc) 



2\k\ 



(66) 



(67) 



One of the key features of this inner product is that 
for an arbitrary wave function the quantity 



(*,«-) = i I da (*;c^* + - 



&|fc|(|*+(fc)| 2 + |*-(fc)| 2 ) , (68) 



is positive definite. 

Once again we shall be interested in calculating the 
probability of the universe in a given time (</>) to have 
a size within the range A or in its complement A. For 
a given initial state ^(cf>,a), we can construct the deco- 
herence functional for a set of histories as proposed in 
Ref. [5] and then take the limit of infinite past 0! — > — oo 
and infinite future 02 — > +oo. 



The off-diagonal terms of decoherence between histories that cross the surface =const. within region A or in A 
is given by 



D(A, A) = / da da' 

J A J A 



(a', 0')§?G (+) (a', 0'; a, 0) *+(a, 0) + (a', 0')^G^ (a', 4>'i <f>)K^- («. 



(69) 



Omitting the negative frequency terms and defining the region A = (— oo, a*) we have 



*;(«', 0')^'GM(«', Ui 0)d * + (a, 0) + 4-+(a, 4>)d^+\a', 0'; a, cj>)d^%{a' , 0') 



D(A, A) = / da / da' 

J —oo J a* 

*;(a', 0')*+(a, cj>)d^G^\a', 0'; a, 0) + G (+) (a', 0'; a, 0)9 ^;(a', 0')^*+(c 
where in our specific case, the Green functions G^' eq. (63) read 

2n\J 2\k\ 



^±ik{v r — v' r ) 



+ 



dk 
2\k\ 



(70) 



(71) 



Let us evaluate Eq. ( 70 ) term by term. The first two terms involving first derivatives of the Green function are zero 
as long as, when 0' = 0, d^G^ (a' , 0'; a, 0) yields a 5{a' — a) and the integrations are in disjoint domains of a. The 
third and fourth terms are more involved and require more attention. Each of them involves eight terms, where four 
of them are null. These are the terms involving either v r or vi only. For instance, the term proportional to 



-co J a* 
oo />a*~ 

/nun 



da / da' / dk' e - ik '<^* + (k') / dk e lkv ^+{k) / dk"\k"\ e lk ' 



dk"\k" 



dv r e l(k ~ k "^ / dv' r e l{k "~ k,) < / dk'** + (k') / dfc*+(/c) 



a* — 



(72) 
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is zero in the limit <f> — > ±00. After the change of variables da — > dv r , it remains a 4> dependence only in the limits 
of integration, hence we can safely take the limit ip — ¥ ±00 before integrating the expression which makes it to go to 
zero. Notwithstanding, the mixed terms that include both vi and v r have a complete different structure. These terms 
are proportional to 



a* />oo 

da I da 1 

-co J a* 



dk' e~ lk '<f*{k') / dk e ikvi ^ + {k)( dk"k" e lk " {a '- a) - dk"k" e tk " {a '- a) 



dk' e- ikW ^* + (k') dke ikv ^+{k)( dk"k" e lfe '>'- Q ) 



—00 




+ / dk'k' e- ik '<^*Ak') / dkk e ikv '^+(k) 



+ I dk'k 1 e- lk ' v '^* + {k') I dkk e lkv ^+{k) 



— OO 
OO 



dk" 

K 
°° dk^_ 

1c" 
K 



a ik (a' — a) 



ik (a' — a) 



dk"k" e k " [a '- a) 

c 

dk" 



k" 



—00 

JW 



dk" 
~kr 



a ik (a — a) 



a ik (a.' — a) 



(73) 



Let us take, for instance, the first term of the first line of the above equation. After some change of variables we 
obtain 



/*oo />oo />oo />oo pO 

/ dk"k" / dvi e -iMk-k") / dfl / e i<(fe"-fe') / dfc'*;(fc') / dk^+(k) e »*( fc - fc ') e <*(*+*') 

/O JO JO Jo J -00 



(74) 



Using again that 



da e lax = ttS(x) ~ ip.v. 



(75) 



one gets four integrals. One of these integral ends up as 

k" 



dk" 



(k" - k')(k" - k) 



(76) 



which has an ultra-violet logarithmic divergence at in- 
finity. Using the same reasoning in the first term of the 
third line, one gets the integral 



dk" 



k"(k" - k'){k" -k) 



(77) 



which now presents an infra-red logarithmic divergence 
the origin. The crucial point is that these are different 
divergencies which cannot cancel each other out. 

In this way, the decoherence functional cannot be made 
diagonal, and hence we cannot construct consistent his- 
tories. 

In fact, one could have anticipated this result. Note 
that the Wheeler-DeWitt equation we are considering is 
completely analogous to the Klein-Gordon equation for 
a massless relativistic particle. However, as pointed out 
in Ref. [S], where the decoherence functional was con- 
structed for a massive relativistic particle, it was observed 
that the off-diagonal terms D(A,A) may become negli- 
gible only if the region A and its complement are much 
larger than the Compton wavelength mT 1 of the particle. 
If we naively take the limit m — ¥ 0, there is no region A 



in which the off-diagonal terms can become negligible. 
Note, however, that the m — > limit of a Klein-Gordon 
particle is tricky and subtle. That is why we have con- 
structed the decoherence functional for the equivalent of 
a massless scalar field from the beginning. 

Note, however, that if one applies the de Broglie-Bohm 
quantum theory to the same problem, one can obtain in- 
formation about the behavior of the early universe, and 
Bohmian bouncing trajectories appear in many circum- 
stances. This was done in detail in Ref. RHZ1. 



VI. CONCLUSION 

We have shown that claims asserting that the Wheeler- 
DeWitt quantization does not eliminate the classical cos- 
mological singularity are not correct without specifying 
the quantum interpretation one is adopting and the quan- 
tization procedure one is taking. 

In fact, there are several papers showing quantum 
bounces in the Wheeler-DeWitt quantization scenario 
that were published much before this sort of claims have 
been first presented (a few examples of the long list of 
papers are [T2lfT6] ). 

This, however, does not diminish in any way the im- 
portance of the series of results obtained in the context of 
loop quantum cosmology. Indeed, loop quantum cosmol- 
ogy has an important advantage over Wheeler-DeWitt 
quantum cosmology inasmuch it has a strong connection 
with loop quantum gravity. Loop quantum gravity has 
much less conceptual problems as a quantum theory of 
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gravity than the canonical quantization procedure that 
leads to the Wheeler-DeWitt equation. Therefore, the 
existence of bouncing solutions in loop quantum cosmol- 
ogy is a significantly relevant result. Note, however, that 
people working in loop or non-loop quantum cosmology 
must be very precise on what quantum theory they are 
taking to interpret their results. 

In connection with the conclusion of Ref . [T] , we have 
shown here that the answers given by the consistent his- 
tories approach are quite fragile. In fact, the existence of 
a quantum bounce strongly depends on the family of his- 
tories one is taking. One can argue that the family with 
only two moments of time, where quantum bounces do 
not exist, encompass the families of histories with three 
moments of time. Take, however, a genuine quantum 
state. In the two-time family we are sure that there is 
no quantum bounce, but in the three-time family this 
question cannot even be posed. This is characteristic 
of the consistent histories approach: the notion of truth 
depends on the family of histories one is taking. This 
ambiguity on the notion of a true statement in the con- 
sistency histories approach can be made quite dramatic 
in other circumstances [3TJ [32] . 

Finally, we would like to stress that the results of this 
paper go much beyond the question about the existence 
of a quantum bounce. It shows that different quantum 
theories may present quite discrepant results when this 
system is the Universe. This finding points out to a 
hope that maybe in cosmology one can find a way do dis- 
criminate between the many proposed quantum theories 
which, asides subjective and philosophical preferences, 
have all the same scientific status in the laboratory. 
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Appendix A: Distributions 



We use throughout the text the following distribu- 
tion(see for instance [33J 131] ) : 



d(x) := I da e tax = n6(x) - ip.v. 



(Al) 



It is the Fourier transform of a Heaviside function, 
pv (-) stands for principal value of -. It solves the fol- 
lowing equation in the sense of distributions 



p.v. [ — ) .;■ = ] , 



(A2) 



The action of the principal value on any function / G 
Co°(R) is given explicitly by: 



lim 



-f(x)dx ■ 



-f(x)dx 



m - f(-x) 



dx. 



(A3) 
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